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Abstract 

A modified de Broglie-Bohm (dBB) approach to quantum mechanics is 
presented. In this new deterministic theory, the problem of zero velocity 
for bound states does not appear. Also this approach is equivalent to stan- 
dard quantum mechanics when averages of dynamical variables are taken, in 
exactly the same way as in the original dBB theory. 
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In the de Broglie-Bohm quantum theory of motion (dBB) which pro- 
vides a deterministic theory of motion of quantum systems, the following 
assumptions are made 0: 

(Al) An individual physical system comprises a wave propagating in space 
and time, together with a point particle which moves continuously under the 
guidance of the wave. 

(A2) The wave is mathematically described by ^f(x, t), a solution to the 
Schrodinger's wave equation. 

(A3) The particle motion is obtained as the solution x(t) to the equation 

. = ft (fJMgjQ m 

with the initial condition x = x(0). An ensemble of possible motions associ- 
ated with the same wave is generated by varying x(0). 

(A4) The probability that a particle in the ensemble lies between the 
points x and x + dx at time t is given by dx. 

In spite of its success as a deterministic theory, this scheme has the draw- 
back that for many bound state problems of interest, in which the time- 
independent part ifj of the wave function is real, the velocity of the particle 
turns out to be zero; i.e., the particle is at rest whereas classically one would 
expect it to move 0. 

To have a better understanding of this problem, let us recall that in the 
the classical Hamilton-Jacobi theory, the trajectory for a particle is obtained 
from the equation 

dS TT ( dS\ 

m +H [ x >irx)= < 2 > 

by first attempting to solve for the Hamilton-Jacobi function S = W — Et 
and then integrating the equation of motion 

dS 

p = rax = — — (3) 

ox 

with respect to time. However, since the classical Hamilton-Jacobi equation 
is not adequate to describe the micro world, one resorts to the Schrodinger 
equation ihd^/dt = —(h 2 /2m)d 2 ^/dx 2 + and solves for the wave func- 
tion The basic idea of dBB is to put \& in the form Re tS ^ h and then 
to identify S as the Hamilton-Jacobi function. The particle velocity is then 
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given by Eq. (|3|) above, and this leads to Eq. ([]]). For wave functions whose 
space part is real, the Hamilton's characteristic function W = constant and 
hence the velocity field is zero everywhere, which is the problem mentioned 
above. 

Another deterministic approach to quantum mechanics, which also claims 
the absence of this problem, is the trajectory representation developed by 
Floyd In one dimension, Floyd uses the generalized Hamilton- Jacobi 
equation 
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where W (' denotes J^) is the momentum conjugate to x. Recently Faraggi 
and Matone [|] have independently generated the same equation (referred 
to as the quantum stationary Hamilton- Jacobi equation) from an equiva- 
lence principle free from any axioms. The general solution to this nonlinear 
equation is given by U 

W = ±(2m) 1/2 (a<f)' 2 + b6 2 + c(f)6)-\ (5) 

where (0, 9) is a set of normalized independent solutions of the associated 
time-independent Schrodinger equation and (a, 6, c) is a set of real coefficients 
such that a, b > 0. Floyd argues that the conjugate momentum W is not the 
mechanical momentum; instead, mx = mdE/dW . The equation of motion 
in the domain [x, t] is rendered by the Jacobi's theorem. For stationarity, 
the equation of motion for the trajectory time t, relative to its constant 
coordinate r, is given as a function of x by |§, f|, [|] 

where the trajectory is a function of (a, 6, c) and r specifies the epoch. Each 
of these non unique trajectories manifest a microstate of the Schrodinger 
wave function, even for bound states. Thus the Bohmian and trajectory rep- 
resentations differ significantly in the use of the equation of motion, though 
they are based on equivalent generalized Hamilton- Jacobi equations. The 
trajectory representation does not claim equivalence with standard quantum 
mechanics in the predictions of all observed phenomena nor does it explain 
the concept of probability density inherent in the Copenhagen interpretation. 
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In this Brief Report, we present a different and modified version of the 
dBB that surmounts the problem mentioned in the introduction. In addition, 
it is closer to the classical Hamilton- Jacobi theory than even the conventional 
dBB. We apply the scheme to a few simple problems and find that it is capable 
of providing a deeper insight into the quantum phenomena. Also it is shown 
that like dBB, the new scheme is equivalent to standard quantum mechanics 
when averages of dynamical variables are taken. 

To this end, we first note that in the dBB, the substitution \I/ = Re lS l h in 
the Schrodinger equation gives rise to two coupled partial differential equa- 
tions, one of which contains an unknown quantum potential term, and this 
leads to a situation apparently different from the classical Hamilton-Jacobi 
theory for individual particles. On the other hand, we note that a change 
of variable \& = Afe lS ^ h in the Schrodinger equation (where M is a constant) 
gives a quantum-mechanical Hamilton-Jacobi equation || 



which is closer to its classical counterpart and which also has the correct 
classical limit. This substitution brings the classical expression (|3|) for the 
conjugate momentum to the form 



We therefore attempt to modify the dBB scheme in a manner similar to the 
latter approach; i.e., we retain assumptions Al, A2 and A4 whereas in A3, 
we identify \l/ = J\fe lS ^ n and use the expression (||) as the equation of motion. 

To see how the scheme works, let us consider the example of the ground 
state solution ty of the Schrodinger equation for the harmonic oscillator in 
one dimension, the space part of which is real. We have 





^ = Af oe - a2x2 / 2 e- lEot/h . 
The velocity field in the new scheme is given by Eq. (j8|) as 



(9) 
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whose solution is 
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X = A ^o?tlm_ 

This is an equation for a circle of radius \ A\ = \x(0) \ in the complex plane 
[Fig. [5(a)]. As usual in such mechanical problems, we take the real part of 
this expression, 

x r = A cos(ha 2 t/m) (12) 

[where it is chosen that at t = 0, x(0) = A is real], as the physical coordinate 
of the particle. It shall be noted that this is the same classical solution for a 
harmonic oscillator of frequency ujq = ha 2 /m. 

We can adopt this procedure to obtain the velocity field, also for higher 
values of the quantum number n. For n — 1, we have 

= A/le^ 2 * 2 / 2 2ax e - lE ' t/h , (13) 

from which 

x = — (~a 2 x + -) . (14) 
im \ x J 

The solution to this equation can be written as 

{ax - l){ax + 1) = Ae 2lha2t/m (15) 

or 

x = + Ae 2iha2t / m . (16) 

Here, the solution is a product of two circles centered about ax = ±1, which 
is plotted in Fig. [j](b) . The physical coordinate of the particle is again given 
by the real part of this expression. For n = 2, the solution can similarly be 
constructed as 

4ax (ax - ^5/2^j * (^ax + ' = Ae 5iha2t/m . (17) 

The trajectory in the complex plane is plotted in Fig. 0(c). Note that in all 
these cases, the velocity fields are not zero everywhere. 

Now let us apply the procedure to some other stationary states which 
have complex wave functions. For a plane wave, we have 
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$ = Ae ikx e -iEt/n^ ^ 
so that x = hk/m and this gives 

x = x(0) + — t, (19) 
m 

the classical solution for a free particle. 

As another example, consider a particle with energy E approaching a 
potential step with height Vq, shown in Fig. ||(a). In region I, we have 



^ = e ikx + Re 



-ikx 



k = J2mE/Ti z 



and in region II, 



1> n = Te***, q = pm{E - V )/h 2 . 
The velocity fields in the two regions are given by 

hk ( e lkx — Re~ lkx \ 



m \ e lkx + Re~ 



ikx 



and 



x n 



hq 



m 



(20) 



(21) 



(22) 



(23) 



respectively. The contours in the complex x-plane in region I, for a typical 
value of the reflection coefficient r = R 2 = 1/2, are given by 



\pl cos 2kxi T — e 



2kx r 



~-e 2kxii = c 



(24) 



[where xi r and xn are, respectively, the real and imaginary parts of xi], and 
are plotted in Fig. @(b). Note that also this case is significantly different 
from the corresponding dBB solution. 

Lastly, we consider a nonstationary wave function, which is a spreading 
wave packet. Here, let the propagation constant k has a Gaussian spectrum 
with a width Ak ~ 1/cr about a mean value k. The wave function is given 

by 



*(x,t) = M 



2na 



1/2 



a 2 + iht/m 



exp 
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2(a 2 + iht/m) 



(25) 
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The velocity field is obtained from (§) as 

% ( x — ia 2 k 



x 



im \a 2 + iht/m J 
Integrating this expression, we get 

. . hk .h x(o) 

x(t) = x(0) + —t + i V 

Separating the real and imaginary parts of this equation (and assuming k is 
real), we get, 



t. 



(26) 



(27) 



x r (t) = x r (0) H 1 H \±t 



m m 



and 



Xi(t) = Xi(0) 



h x r (o) 



m a* 



t. 



(28) 



(29) 



For the particle with x(0) = 0, we obtain x r {t) = (hk/m)t and Xj(t) = 0; 
i.e., this particle remains at the center of the wave packet. Other particles 
assume different positions x(t) at time t as given by the above expression, 
which indicates the spread of the wave packet. 

A new element in the present approach to quantum mechanics (though 
it is quite familiar in elementary mechanical problems) is that the particle 
coordinates are having real and imaginary parts. Similar is the case for the 
velocity, given by Eq. (||). If we are interested in the real part of this velocity, 
then one can write 



dx \ dx ' 



Xr ~2im V*V ' [6()) 

an expression quite analogous to the expression for velocity field in the dBB 
scheme (though not identical with it). One can also write 



dx r 



• dx r 



2im 

where use is made of the fact that for complex derivatives, ^ = 
equation is identical to the dBB velocity field at all points x = x r . 



(31) 



This 
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Now consider an ensemble of particles, whose initial density function along 
x = x r is given by 

P(x r , t = 0) = ^*(x = x r ,t = 0)tf (a; =x r ,t = 0). (32) 

Now let us ask whether P(x r ,t) = ty*(x = x r ,t)^(x = x r ,t) be identified as 
the density function along x = x r at any time £, given the velocity field (||). 
For this to be in the affirmative, the conservation equation 

dP d 

m + a^> = (33) 

must hold, with x r given by Eq. ( |3~0| ) above. It is easy to see that this is 
true, provided V(x = x r ) is real. Thus one can retain the assumption (A4), 
with ty*(x = x r ,t)ty(x = x r ,t) as the density function at any time t. This 
guarantees that the averages of dynamical variables O constructed using 
the real variables x r and t, computed over the measure P will necessarily 
agree with the quantum mechanical expectation values of the corresponding 
hermitian operators O at all future times t; i.e., 

0=1 POdx r = / ^*(x = x r ,t)0^(x = x r ,t)dx r , (34) 

J t=constant J 

in the same way as in dBB. Thus the new scheme is equivalent to standard 
quantum mechanics when averages of dynamical variables are taken, just as 
in the case of the original dBB approach. 

Generalization of the formalism to more than one dimension is not at- 
tempted in this Brief Report. Its application to those other physical problems 
of interest and a host of other issues, which need careful analysis, shall be 
addressed in future publications. 

Finally, we note that the dBB identification \l> = Re lS ^ n does not help to 
utilize all the information contained in the wave function while solving the 
equation of motion mx = dS/dx, though it provides a conservation equation 
for R 2 . In the present formulation, we obtain a quantum Hamilton- Jacobi 
equation, which is closer to the classical one and which does not have any 
exotic quantum potential term. But still we could get the conservation equa- 
tion for probability density, as demonstrated above. The positive results we 
obtained for the harmonic oscillator and potential step problems themselves 
are indicative of the deep insight obtainable in such problems by the use of 
the new scheme. The price we had to pay for this is the appearance of an 
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imaginary part in the position and velocity of the particle, but the technique 
is considered to be quite normal in elementary mechanical problems (and 
not yet in quantum mechanical ones). Therefore, it is desirable to further 
explore the consequences of the new scheme. 
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Figure 1: Trajectories in the complex X-plane for the harmonic oscillator 
(X = ax), (a) The n = case, where contours are plotted for X(0) =1, 2, 3 
and 4. (b) The n = 1 case, where contours are plotted for X(0) =1.2, 1.35, 
1.45 and 1.55. (c) The n = 2 case, where contours are plotted for X(0)=1.8, 
1.9, 2.0 and 2.1. 
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Figure 2: (a) Potential step, with V = for x r < 0, V = Vq for x r > 
and energy E > Vq. (b) Trajectories in the complex plane for particles 
approaching the potential step. Contours for c = -4, -3, -2, -1 and are 
plotted. 



